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ABSTRACT 

We present an empirical connection between cold gas in galactic halos and star forma- 
tion. Using a sample of more than 8,500 Mgii absorbers from SDSS quasar spectra, 
we report the detection of a 15 ct correlation between the rest equivalent width Wq 
of Mg II absorbers and the associated [O ii] luminosity, an estimator of star formation 
rate. 

This correlation has interesting implications: using only observable quantities we 
show that Mg ii absorbers trace a substantial fraction of the global [O ii] luminosity 
density and recover the overall star formation history of the Universe derived from 
classical emission estimators up to z ~ 2. We then show that the distribution function 
of Mgii rest equivalent widths, diV/dWo inherits both its shape and amplitude from 
the [On] luminosity function These distributions can be naturally connected, 

without any free parameter. 

Our results imply a high covering factor of cold gas around star forming galaxies: 
C > 0.5, favoring outflows as the mechanism responsible for Mgii absorption. We 
then argue that intervening Mg ii absorbers and blue-shifted Mg ii absorption seen in 
the spectra of star forming galaxies are essentially the same systems, implying that 
the observed outflowing gas can reach radii of ~ 50 kpc. These results not only shed 
light on the nature of Mg ii absorbers but also provide us with a new probe of star 
formation, in absorption, i.e. in a way which does not suffer from dust extinction and 
with a redshift-indepcndent sensitivity. As shown in this analysis, such a tool can 
be applied in a noise-dominated regime, i.e. using a dataset for which emission lines 
are not detected in individual objects. This is of particular interest for high redshift 
studies. 

Key words: absorbers: Mgii ~ star formation rate - quasars - outflows 



1 INTRODUCTION 

In the past twenty years, we have seen enormous advances in 
our understanding of galaxy evolution. Gravitational lensing 
has allowed us to interpret the properties of the host dark 
matter halos, population synthesis models can characterize 
the emission properties of galaxies over 3 orders of magni- 
tude in stellar mass, and the star formation history of the 
Universe has been explored up to high redshift, revealing a 
peak at z ~ 2. Along with this, the importance of feedback 
processes has emerged and one of the central questions which 
now needs to be addressed is how galaxies accrete, process 
and return gas into the inter-galactic medium (IGM). 



Absorption line spectroscopy provides us with a pow- 
erful tool to address this question, to explore in detail the 
distribution of gas around galaxies and to study its inter- 
play with star formation. It gives us access to a wide array 
of elements and allows us to detect low-density gas that is 
orders of magnitude below the detection threshold of most 
other techniques. 

Historically, the study of QSO absorber-galaxy relation- 
ships began with Mg ll absorption systems. This choice was 
mainly driven by practical constraints: among the dominant 
ions in H I gas only Mg ll could be studied from the ground 
at low redshifts {z < 2). Following the first confirmation of 
a Mgll absorber-galaxy connection by Bergeron (1986), a 
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considerable amount of work has focused on revealing the 
nature and origin of these systems. Mgii is known to be a 
tracer of galactic halos, routinely seen up to about 50kpc 
around galaxies. However, despite significant observational 
efforts in the last two decades, fundamental questions re- 
garding the physical nature of the absorbing gas remain 
unanswered. Various models have been proposed in the lit- 
erature: outflows (Bond et al. 2001), infalling gas (Tinker 
& Chen 2008), extended disks (Steidel et al. 2002), etc. but 
no consensus has been reached. Observationally, no com- 
pelling correlations between absorber properties and those 
of their associated galaxies have been reported. In particu- 
lar, the physical process(es) governing the observed equiva- 
lent widths Wo, the basic observable quantity, and the origin 
of its distribution function are not yet understood. 

In this paper we present the detection of a correla- 
tion between the observed rest equivalent width of Mgii 
absorbers and the associated [O ll] luminosity and show that 
this relation has important implications. First, Mgll gas ap- 
pears to trace a substantial fraction of the [O ll] luminosity 
density of the Universe. These absorber systems can there- 
fore be used as a new probe for star formation over cosmolog- 
ical times. We show that the observed redshift dependence 
reproduces the decline in cosmic star formation rate since 
z ~ 2 observed with emission line surveys. Finally, we de- 
rive a simple relation between the probability distribution 
function (p.d.f.) of Wq and Loii and show that these two 
quantities are related. 

This study was motivated by a number of recent re- 
sults, indicating a possible link between the presence of 
strong Mgll absorbers and the star formation rate of their 
associated galaxies: Zibetti et al. (2007) analyzed the col- 
ors of stacked images of quasars with strong Mg ii absorbers 
and reported that the best-fit spectral energy distributions 
(SED) require stronger emission lines for stronger Mgll ab- 
sorbers. Wild et al. (2007) created composite spectra of a 
subset of Mgii-selected absorbers and was able to detect 
their mean [O ii] emission. Rubin et al. (2009) studied the 
properties of a Mg ll absorber-galaxy association and showed 
that the absorbing gas found at about 15 kpc from the galaxy 
is consistent with having been ejected during the last burst 
of star formation of this galaxy. Finally, Nestor et al. (2010, 
in prep) studied the properties of two ultra-strong Mgll ab- 
sorbers {W{) = 5.6 and 3.6A) and showed that the associ- 
ated galaxies are in the starburst and post-starburst phase, 
respectively. 

Using a different approach, by studying self-absorption 
in galaxy spectra, Weiner et al. (2009) showed that all types 
of star forming galaxies trigger outflows seen through blue- 
shifted Mgll absorption. Such an analysis does not provide 
any information on the distance between the absorbing gas 
and the galaxy, and the connection to intervening absorber 
systems is not straightforward. However, we will show that 
our new results suggest a strong link between these two trac- 
ers of gas in and around galaxies. 

The outline of this paper is as follows: after briefly re- 
viewing useful properties of Mg ii absorbers and [O ii] emis- 
sion as a tracer of star formation, we describe the spectro- 
scopic data analysis in Section 3. We present the correlation 
between [O n] luminosity surface density and Mg ii equiva- 



lent width in Section 4. Section 5 presents the connection 
between 3>(Loii) and the distribution of Mgii rest equiv- 
alent widths. Finally, we discuss the implications of these 
results in Section6 and summarize our work in Section 7. 

Throughout this paper we use the cosmology: Q,m = 
0.3, Oa = 0.7 and Ho = 70km/s/Mpc. 



2 ASTROPHYSICAL PRELIMINARIES 

2.1 Properties of Mgii absorbers 

Strong Mgll absorbers, classically defined with Wo'^'''^^ > 
0.3 A, are known to be associated with a range of galax- 
ies, with 0.1 < L < 10 L* (Bergeron & Boisse 1991; Stei- 
del, Dickinson, & Persson 1994; Nestor et al. 2007) and are 
found at impact parameters ranging from a few to more 
than 50 kpc (Zibetti et al. 2007; Steidel et al. 1997; Tinker 
& Chen 2008). In the range Wq^''^'^ > lA Mgll absorp- 
tion lines are in general saturated (Nestor et al. 2005) and 
no column density information can be directly extracted 
from the observed equivalent width, which can instead be 
used as an estimate of the overall gas velocity dispersion. 
High-resolution spectroscopic observations reveal that the 
absorption often originates from several velocity components 
(Churchill, Kacprzak, & Steidel 2005). Empirically, it has 
been found that in this regime the velocity dispersion of the 
gas follows Av ~ 120(Wo/A)km/s (see Ellison 2006, Fig. 
3). 

A fundamental quantity describing the statistical prop- 
erties of these absorbers is the distribution of Mgii rest 
equivalent width per unit redshift, d^N/dWo dz. This quan- 
tity has been accurately characterized using the Sloan Digi- 
tal Sky Survey (SDSS) by Nestor et al. (2005). These authors 
found that it is well described by an exponential distribu- 
tion: 

dWo W* ^ ' 

with the maximum likelihood values N* = 1.00 ± 0.13 (1 -|- 
^-)0.22±o.i7 ^* ^ 0.44 ±0.03 (l + z)°-83±o io. This result 
is valid down to a rest equivalent width limit of ~ 0.3 A, 
below which the distribution follows a power-law behavior, 
making the overall distribution well described by a Scliechter 
function. So far, no compelling model has been able to ex- 
plain the origin of this distribution. 



2.2 [O ii] emission as a tracer of star formation 

Over the past 20 years many techniques have been explored 
to estimate the star formation rates (SFRs) of galaxies and, 
all together, have established the steep rise in the SFR den- 
sity of the Universe from the present epoch to z ~ 2 (Hop- 
kins & Beacom 2006 and references therein). One of the best- 
understood SFR indicators is Ha A6563, whose luminosity 
is directly proportional to the hydrogen-ionizing radiation 
from massive stars (M> 10 M©) and therefore provides a 
near-instantaneous (r < 10 Myr) measure of the SFR with 



© 0000 RAS, MNRAS 000, 000-000 



Probing star formation with absorption line systems 3 




1000 



100 



10 



n I I I r 




1 I . . I 



_l I I L_ 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
redshift 



0.7 1 2 3 4 5 6 

wr^(Mgll) [ A ] 



Figure 1. Left: redshift distributions of detected Mgll absorbers (red) and background quasars with z < 3 (gray). The filled histogram 
represents the redshift range for which both Mg II and [O ll] lines can be detected. Right: observed distribution of MgII-A2796 A rest 
equivalent widths. Below 1 A the completeness becomes significantly lower than unity and N(VFo) departs from an exponential distribution. 



minimal dependence on the physical conditions of the ion- 
ized gas (Kennicutt 1998). Such observations are however 
limited to z < 0.4 in the optical. At higher redshifts, the 
[O ll] A3727 nebular emission line has been used as an alter- 
native SFR indicator (e.g. Gallagher, Hunter, & Bushouse 
1989; Kennicutt 1992; Hopkins et al. 2003; Kewley et al. 
2004; Mouhcine et al. 2005). 

The [O ll] emission doublet is the shortest-wavelength 
strong emission line from low-density, photoionized galactic 
nebulae. The nebular lines eflfectively re-emit the integrated 
stellar luminosity short- ward of the Lyman limit. Only stars 
with masses M > 10 Mq and lifetimes r < 20 Myr con- 
tribute significantly to the integrated ionizing flux, so the 
emission lines provide a nearly instantaneous measure of the 
SFR, independent of the previous star formation history. 
[O ll] luminosity depends explicitly on the chemical abun- 
dance and excitation state of the ionized gas, and suffers a 
larger amount of dust extinction than Ha. Therefore, unlike 
the Balmer recombination lines, [O ll] is not directly pro- 
portional to the SFR and must be calibrated either empir- 
ically (Kennicutt 1992; Kewley et al. 2004) or theoretically 
(Barbaro & Poggianti 1997; Chariot & Longhetti 2001). On 
average, it has been found that 



SFR [Mo/yr] ~Ax 10"*' Loii [erg/s] 



(2) 



The value of A is of order unity, with the precise value de- 
pending on the broad band luminosity and dust content of 
the galaxy, etc. (Moustakas, Kennicutt, & Tremonti 2006). 



3 DATA AND ANALYSIS 

Strong Mg ll absorbers can be detected in SDSS quasar spec- 
tra over the redshift range 0.4 < z < 2.2. The spectroscopic 



observations make use of 3" fibers receiving photons origi- 
nating from the entire path to the quasars (convolved with 
the seeing of the observations). The detection of a strong 
Mgii absorber indicates the presence of a galaxy close to 
the line-of-sight. If this galaxy is starforming, its narrow 
emission lines may be identified on top of the quasar spec- 
tral energy distribution. While the direct detection of these 
lines is possible only for a small number of absorber systems 
at low redshift (see Noterdaeme, Srianand, & Mohan 2009), 
a statistical approach enables the measurement of the mean 
luminosity of such lines well below the noise level of individ- 
ual spectra and allow us to study the relationship between 
emission and absorption lines. 



3.1 The dataset 

The analysis makes use of a sample of Mgll systems com- 
piled by Quider et al. (2010) using the method presented in 
Nestor et al. (2005) and extended to the SDSS DR4 dataset 
(Adelman-McCarthy et al. 2006). In this section we briefly 
summarize the main steps involved in the absorption line 
detection procedure. We refer the reader to Nestor et al. 
(2005) for more detail. 

For each quasar spectrum in the SDSS DR4 database, 
Quider et al. (2010) fit a pseudo-continua in an iter- 
ative fashion using a combination of cubic-splines and 
Gaussians for both emission and absorption features. The 
continuum-normalized SDSS QSO spectra were searched 
for Mg IIAA2796, 2803 doublets using an optimal extraction 
method employing a Gaussian line-profile to measure each 
rest equivalent width Wq. All candidates are interactively 
checked for false detections, a satisfactory continuum fit, 
blends with absorption lines from other systems, and spe- 
cial cases. The identification of Mg 11 doublets required the 
detection of the 2796 line and at least one additional line, 
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Figure 2. Median flux in erg/s/cm^/A from composite spectra 
in the region of the [O ll] emission line, for five bins of Mg II rest 
equivalent width. The red line is the best fit emission doublet. 
The spectra have been offset for display purposes. 



fraction of direct detections can be made at low redshift 
and will be discussed in future work. In the present anal- 
ysis we estimate the average [O ll] emission of samples of 
Mgll absorbers by creating composite spectra in the ab- 
sorber rest frame. In this procedure masks are applied to 
absorption lines expected at the redshift of all known Mgll 
absorber systems as well as sky emission lines. The contin- 
uum is defined with an iterative running median of sizes 
ranging from 500 to 15 pixels. Narrow features caused by 
unidentified systems with z 7^ Zabs cannot be masked with 
a reliable uniformity, and are thus treated as an additional 
noise contribution. Each continuum-subtracted spectrum is 
converted into units of luminosity surface density (see be- 
low) and samples of spectra are combined to create average 
and median composites. In general the results from these 
two estimators are found to be similar. Since the former 
ones are usually noisier we will focus only on median values 
below. 

In each composite spectrum the [O ll]AA3727,3730 dou- 
blet was fit with a double Gaussian with line ratio fixed at 
4:5 and freely varying line width. Line luminosity surface 
densities (see below) were measured from the fitted Gaus- 
sian parameters. In order to estimate the noise level of the 
line detection we applied the same fitting procedure to the 
region of the spectrum 3500 < \reat < 3600 A where no 
strong emission/absorption line is expected. 

Examples of average composite spectra around the [O 11] 
line are shown in Fig. 2 for samples of Mg 11 absorbers with 
different rest equivalent widths. We note that no shift in the 
centroid of the [0 11] lines with respect to those of the Mg 11 
lines is detected. These results indicate a strong correlation 
between [0 11] emission and Mg 11 absorption rest equivalent 
width. This correlation will be analyzed quantitatively be- 
low. 



the most convenient being the 2803 doublet partner. A 5a 
significance level was required for all A2796 lines, as well as a 
3(7 significance level for the corresponding A2803 line. Only 
systems Q.lc blue- ward of the quasar redshift and red- ward 
of Ly— a emission were accepted. From this catalogue we 
select systems where both Mgii and [O ii]AA3727,3730 are 
simultaneously covered by the accessible wavelength range. 
The final sample contains 8523 absorption line systems with 
0.7 <Wo< 6.0 A and OA < z < 1.3. 

The redshift distribution of the systems used in this 
analysis is presented in the left panel of Figure 1 as the filled 
histogram. The right panel shows the observed distribution 
of absorber rest equivalent width Wq. As described above, 
the intrinsic distribution of rest equivalent widths follows an 
exponential distribution down to Wo — 0.3 A (see Eq. 1). 



3.2.2 Aperture correction 

The available sample of absorbers for which both Mg 11 and 
[O 11] lines are simultaneously accessible spans a substantial 
redshift interval: from 0.4 to 1.3. In this range, the physical 
area corresponding to a given angular aperture varies by a 
factor of about 2.5. This is illustrated in Fig. 7 in the Ap- 
pendix in which we show the physical size covered by the 
3" diameter SDSS fibre aperture as a function of redshift. 
In order to combine coherently, or compare, signals from 
different redshifts, we estimate the [0 11] luminosity surface 
density, Slqh- For each Mgll absorber, we thus convert the 
spectrum from fiux units into luminosity units, at the red- 
shift of the absorber. We then divide by the surface area of 
the SDSS fibre at this redshift, to obtain luminosity surface 
density: 



3.2 Measuring [O 11] luminosities 

3.2.1 Composite spectra and line measurement 

In general, the [O 11] emission lines from galaxies with z > 
0.4 are not detected in individual SDSS spectra. A small 



Sloii(z) 



Loujz) 



(3) 



where — vr^y and Of is the angular radius of the SDSS 
fibre and Da [z) is the angular diameter distance at the red- 
shift of the absorber. 
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Figure 3. Left: mean [On] luminosity surface density, Slqjj (averaged within the 3-arcsec fibre) as a function of Mgll rest equivalent 
width Wo. The solid line is the best fit power-law. Right: mean [On] luminosity surface density measured in three redshift bins, as a 
function of Wq. The dashed lines show the best 3-parameter fit to the entire sample. 



Mgii ABSORPTION AS A TRACER OF 
STAR FORMATION 



provide us with a powerful probe of star formation over cos- 
mological times. This will be discussed in more detail below. 



In this section we use Mgll absorbers as tracers and show 
they can be used as a probe of star formation. This is done 
independently of any relation between Mgll absorbers and 
galaxies. 



4.1 Mgll - [Oil] correlation 

Using a logarithmic binning in Mgll rest equivalent width, 
W^^^^, we create composite spectra in the absorber rest- 
frame and estimate Slqh , the median [O ll] luminosity sur- 
face density. The observed relation between these two pa- 
rameters is presented in Fig. 3 with red data points. Mgll 
absorbers with equivalent widths 0.7 < Wo < 6 A are seen to 
strongly correlate with [O ll] emission, spanning more than 
a factor 30 in luminosity surface density. The observed re- 
lation, averaged over the redshift range 0.4 < z < 1.3, is 
simply described by 



(Elo„(W^o)) 



Wo 

lA 



(4) 



with A = (1.48 ±0.18) x 10^'^ ergs"^ kpc^ and a = 1.75 ± 
0.11. This fit is shown in Figure 3 with the solid line. 

This result shows that, without any prior knowledge of 
the underlying galaxy distribution, Mgll absorbers can be 
used to trace the distribution of [On] emission, i.e. star for- 
mation. This property is interesting since the detectability 
of absorber systems does not depend on redshift and is not 
significantly affected by dust extinction. They can therefore 



4.1.1 Redshift dependence 

In order to investigate the redshift dependence of the above 
relation, we measure the mean [0 11] luminosity surface den- 
sity for three redshift intervals: 0.36 < zi < 0.67 < 22 < 
0.99 < Z3 < 1.3. The results are shown in the right panel of 
Fig. 3 where each color corresponds to a given range of Mg 11 
rest equivalent widths, resulting in 12 independent data 
points. Remarkably, they appear to be distributed regularly 
in the Loii — z plane. This property further demonstrates 
that the relation between Mgll absorption properties and 
observed [0 11] luminosity surface density holds over a sig- 
nificant cosmological time span. This figure illustrates that 
Mgll absorbers are found, on average, in regions having a 
predictable value of [0 11] luminosity surface density. This 
correspondence is valid over a large range in both redshift 
and [0 11] luminosity. 

A global fit to all the data points in the redshift-Wo 
plane gives: 



(Slo„(W^o,2)>=A 



(i + zY 



(5) 



with A = (7.9 ± 2.0) x 10^® ergs"^ kpc^ a = 1.7 ± 0.1 and 
/3 = 1.1 ± 0.4. The results of this global fit are shown in 
the right panel of Fig. 3 as dotted lines. Using this simple 
formula, we are able to fit 12 data points with 3 parameters 
resulting in a reduced of 1.3. The redshift dependence of 
this relation is remarkably weak. The mean [O 11] luminos- 
ity surface density changes by less than 50% over a unity 
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Figure 4. Left: Relative contribution of [On] luminosity surface density, Slqjj as a function of Mgn rest equivalent width Wq. The 
solid line is the best fit given by Eq. 9. Right: the [On] luminosity volume density from Mgn absorbers (with Wo > 0.7 A in orange and 
all systems in red) computed at three different redshifts. The solid line shows the best fit trend computed within the range for which 
the incidence of Mgll absorbers (dTV/dz) is measured. The other data points show direct [On] luminosity density measurements from 
the literature. Note that no dust correction has been applied to any of these [O n] luminosity measurements, allowing direct comparison 
between samples. Mg n absorbers appear to trace a substantial fraction of the global [O n] luminosity density in the Universe. 



redshift range, corresponding to a timescale of about 5 Gyr. 
Understanding the origin of this dependence would be of 
particular interest and would allow us to directly use Mgll 
absorbers as a probe of star formation without the need of 
detecting [O ll] emission. 



(with and without absorbers). In the second expression, 
D\{z) (1 + z)'^ dO, is the co-moving surface of the aperture, 
AN /Ar is the mean number of absorbers expected within 
a pathlength dr and AL is the [O ll] luminosity observed 
within the corresponding volume element. Introducing the 
cosmology-dependent function 



4.2 Oil luminosity density 

The SDSS fibers used to observe quasars can be considered 
as a set of cones integrating light emitted along the line-of- 
sight. Each quasar spectrum therefore probes a given volume 
and the set of (high-redshift) quasar spectra can be con- 
sidered as a random sampling of the low-redshift Universe. 
Within this volume, we estimate the [O ll] luminosity den- 
sity. Con, traced by Mgn absorbers. Note that we are only 
interested in estimating a density and therefore do not need 
to include any emission contribution originating outside the 
fibers. We then use Con to estimate the star formation rate 
density probed by these systems as a function of redshift 
and compare it to the total value in the Universe. 

The co-moving [O ll] luminosity density probed by Mg ll 
absorbers is given by 

1 AL AN 



Coii{z) 



{1 + zY D\{z)AQ. dr 



(6) 



where L is the [O ll] luminosity associated with absorbers 
observed within the co- moving volume Vc, which is accessi- 
ble along the pathlength towards the background quasars 



£{z) 



Az 



dr {1 + zY 



Ho V!^M(l + ^)^ + nA 

c {l + zY 



(7) 



we can simply express the luminosity density probed by a 
population of Mgll absorbers as a function of observable 
quantities 



Coii(z) = 



(8) 



where AN/Az is the line-of-sight incidence of absorber sys- 
tems and Sloii the mean luminosity surface density associ- 
ated with them. Eq. 8 indicates that, for a given cosmology, 
the [On] luminosity density traced by Mgll absorbers can be 
derived without any assumption. In particular it does not 
depend on the spatial extent of the gas and the relation 
between absorbers and galaxies. A similar formalism was 
used by Wolfe, Gawiser, & Prochaska (2003) and Wild et 
al. (2007) to estimate the contribution of damped Lyman- 
a systems and Call absorption systems to the overall star 
formation rate. 

In order to account for the varying completeness of the 
absorber sample with Wo, we estimate Eq. 8 by integrating 
over the intrinsic rest equivalent width distribution 
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Coii{z) = £{z) / dWo 



dWodz 



SLon(W^0,2) 



(9) 



where d^N/dWodz is the completeness-corrected distribu- 
tion of rest equivalent widths given by Nestor et al. (2005), 
see Eq. 1. The relative contribution of the integrand is shown 
in the left panel of Fig. 4 for the sample spanning the entire 
redshift range. It shows that most of the [O ll] luminosity 
surface density, and therefore most of the star formation, is 
associated with systems having Wo ~ 2 A. Interestingly, the 
same result was obtained regarding the dust content of Mg II 
absorbers (Menard et al. 2008). We will return to this point 
in the discussion. We note that estimating the integral in 
Eq. 9 from the lowest rest equivalent width available in this 
analysis, i.e. Wo ~ 0.7 A, gives a result only 15% lower than 
integrating from zero with the extrapolated fitting functions 
at the mean redshift of the sample. 

We compute Lou using Eq. 9 for the three redshift 
bins defined above and the results are given in Table 1. The 
results are presented in the right panel of Figure 4. The or- 
ange data points show direct measurements for our sample 
of Mgll absorbers with Wo > 0.7 A and the red ones show 
the estimate for all Mgll absorbers, i.e. including the con- 
tribution from weaker systems. For comparison, we show a 
compilation of measurements of the total [O ll] luminosity 
density obtained from narrow band filter and emission line 
surveys. Note that no dust correction has been applied to 
any of these datasets, allowing direct comparisons between 
observed quantities. This figure shows that, at z ~ 1, the 
amount of [On] luminosity density traced by Mgll absorbers 
is a substantial fraction of the total density in the Universe 
estimated from direct emission measurements. Given the 
non-negligible scatter in the reported values of Con from 
photometric and spectroscopic surveys at around 2: = 1, it 
is not yet possible to precisely quantify the fraction of [O 11] 
luminosity traced by Mgll absorbers. Nevertheless our re- 
sults show that at least half of the [O 11] luminosity density 
is traced by strong Mg 11 absorbers. This implies that a high 
fraction of star-forming galaxies are associated with strong 
Mgll absorbers and vice- versa. In other words, the presence 
of Mg II absorbers around a galaxy depends on its (recent) 
star formation rate. 

To further illustrate this connection, we explore the red- 
shift dependence of Con- As shown above, this quantity is 
given by the product of two observables: d'^N /dWo dz, which 
has been accurately measured up to z = 2.2 (see Eq. 1) and 
Eloii(Wd,2) which is parametrized in Eq. 5 up to 2; = 1.3. 
In order to estimate Con up to 2 = 2.2, we use an extrapola- 
tion of Eq. 5, a relation which is weakly redshift dependent. 
The result is shown in Fig. 4 with the solid line. The over- 
all redshift dependence appears to be in agreement with the 
trend given by direct [O 11] emission surveys. 



4.3 Star formation history 

We now show how the global star formation probed by Mg 11 
systems compares to the overall star formation history in 
the Universe estimated by various techniques. As mentioned 
above, [O 11] luminosity can be used as an estimator of star 



redshift 


Con for VKo > O.7A 


Con for all Wo 


interval 


[1039 crgs-i Mpc-3] 


[1039 ergs-l Mpc-3] 


0.36 < z < 0.67 


1.08 ±0.13 


1.37 ±0.16 


0.67 <z< 0.99 


1.35 ±0.17 


1.64 ±0.21 


0.99 < 2 < 1.30 


1.70 ±0.24 


2.00 ±0.28 



Table 1. [On] luminosity density traced by Mgll absorbers as a 
function of redshift. No dust corrections have been applied. 
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Figure 5. Star formation rate density as a function of redshift. 
The gray data points are from the Hopkins (2004) compilation 
and use various types of estimators. The orange and red data 
points show the SFR density traced by Wo > 0.7A and all Mgll 
absorbers. The red curve shows the SFR density estimate using 
dN/dz measured up to 2 = 2.2 and the extrapolated [Oll]-Mgll 
relation. 



formation rate, however the scaling between these two quan- 
tities is not straightforward: it is based on empirical rela- 
tions, calibrated at various redshifts and subject to dust ex- 
tinction corrections (Kennicutt 1998; Moustakas, Kennicutt, 
& Tremonti 2006). Here we choose to apply the average scal- 
ing coefficient used in a recent analysis by Zhu et al. (2009) 
for galaxies with 0.75 < z < 1.45 (see their Table 2). These 
authors made use of the empirical correlation derived by 
Moustakas et al. (2006) between the absolute B-band mag- 
nitude, and the L[Oll]/SFR ratio. This calibration statisti- 
cally accounts for the gross systematic effects of reddening, 
metallicity, and excitation, all of which correlate with opti- 
cal luminosity, and has been shown to work reasonably well 
for star-forming galaxies at 0.7 < z < 1.4 (Moustakas et al. 
2006; Cooper et al.2008). Averaging their scaling coefficients 
over the redshift interval 0.75 < z < 1.45, we get 



= 4.2 10"''^ X L 



on M0 yr" 



Mpc" 



(10) 
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where Loii is in erg s~^. With the hmitations of such a 
scaling in mind, we convert our [O ii] luminosity density es- 
timate into a star formation rate density and present the 
results in Figure 5. The gray data points show the estimates 
summarized in the Hopkins (2004) compilation'^. We show 
the estimated star formation rate density from Mgii ab- 
sorbers with Wq > 0.7 A and all systems with the orange 
and red data points respectively. The solid line shows the 
best fit trend obtained in Eq. 9, using the [O ii] - Mg ii scal- 
ing given in Eq. 5 and dN/dz measurements available up to 
z = 2.2. As can be seen, both the overall amplitude and the 
redshift dependence arc in good agreement with other star 
formation rate estimators. The global star formation probed 
by Mgll absorbers can therefore be used as an independent 
probe of the star formation history in the Universe. 

Interestingly, extrapolating the comparison between the 
SFR density traced by Mgll absorbers and observations of 
the global SFR density to higher redshifts suggests that the 
product Sloii ^ dN/dz (see Eq. 8) should reach a max- 
imum value at around redshift two and then decrease at 
higher z. Since Sloii(2, W5) appears to be weakly redshift 
dependent, the expected decrease in this product should be 
driven by a decrease in dN/dz beyond redshift two. This can 
be tested observationally using near-infrared spectroscopy of 
an ensemble of quasars, with instruments such as X-shooter 
(D'Odorico et al. 2006) or the upcoming FIRE spectrograph 
Simcoe et al. (2008). 



5 THE [O II] LUMINOSITY FUNCTION AND 
THE DISTRIBUTION OF Mgii REST 
EQUIVALENT WIDTHS 

In the previous section we made use of Mgll absorbers as 
tracers and showed that, independently of their relations to 
galaxies, they can be used as a probe of star formation in 
the Universe. The connection between Mgii absorbers and 
star formation can be explored one step further. We now 
show that the distribution of Mgll rest equivalent widths 
dA''/dWo, and the [On] luminosity function, i.e. a measure 
of the jjrobability distribution function (p.d.f.) of star forma- 
tion rate in the Universe, are connected, both in amplitude 
and shape. Our goal here is only to reveal the existence of 
this relation and explore its implications. A more detailed 
comparison and study of possible deviations is left for future 
work. Moreover, as explained below, the scatter in the re- 
ported values of the [0 11] luminosity function (see Fig. 5 and 
6) currently limits the level to which we can meaningfully 
refer to the amplitude of the [0 11] luminosity function. 

Given the existence of a relation between the mean [O 11] 
luminosity surface density and Wq we can express, in a sta- 
tistical sense, the [0 11] luminosity function $(Loii) in terms 
of Wo (Mgll). We can write: 



d^iV dWo £{z) 



dWndz dlogL Ca 



(11) 



$(logL) = 



d^N 



dlogLdV C 



1 The more recent compilation by Hopkins et al. (2006) focuses on 
star formation measurements derived from rest-frame UV obser- 
vations rather than emission-line based star formation estimates. 
It is therefore less appropriate for our comparison. 



where L represents the mean [0 11] luminosity per galaxy 
traced by a strong Mg 11 absorber, a is the cross-section for 
absorption and C is the covering factor for Mgll absorp- 
tion within this cross-section, which can be expressed as the 
product of the fraction of galaxies giving rise to Mgll ab- 
sorption and the spatial cross-section taking into account the 
angular variations and dumpiness of the gas: C = CfCn- 
The function £{z) carries the dependence on cosmology, as 
defined in Eq. 7. 

In order to show the connection between the distribu- 
tion of Mg II rest equivalent widths and the [O 11] luminos- 
ity function, we use a simplifying cissumption: we consider 
that the cold gas distribution around star forming galaxies 
is isothermal, as motivated by Chclouche et al. (2008). In 
this case, the 2-dimensional density distribution goes like 
p2D oc and the probability of observing a system at an 
impact parameter 6, p2Y){b)bdh, is constant up to a maxi- 
mum impact parameter ftmax- The uniform cross-section is 
then simply given hy o — n femax- The fraction of galaxy light 
falling within the fiber of projected radius rf is, on average, 
given by (r//i)niax)^ and the mean [0 11] luminosity per Mgll 
galaxy is Loii,gai = Slqh f, where Slqh is the mean [0 11] 
luminosity surface density derived above. This allows us to 
express the luminosity function as: 



$(logL) = 



d^N dWo £{z) 
dWodz dlog(ELo„) Co- 

d^JV £{z) 
dWodz a Ca ' 



(12) 



with a = dlogELoii/dW^o- Using this approximation, we 
now compute the [0 11] luminosity function expected from 

Mgll absorbers, as a function of Wq. The first two terms of 
the above equation are observable quantities, described by 
Eq. 1 and 5. Various observations of Mgii absorbers-galaxy 
associations indicate that the gas extends up to projected 
radu of about 50 kpc (Stcidel et al. 1995, Zibetti et al. 2007, 
Chen & Tinker 2008). We wiU therefore use &max = 50 kpc 
to characterize the cross-section in our estimate of Eq. 12. 

We compute Eq. 12 at z — 1.2, a redshift at which 
several direct measurements of the [On] luminosity func- 
tion are available (see below). We first use C = 1, implying 
that all star forming galaxies are surrounded by Mgll gas. 
The resulting [0 11] luminosity function is shown in the left 
panel of Figure 6 with the solid red line. For comparison, 
we show several estimates of the [0 11] luminosity function 
at the same redshift from narrow-band imaging surveys of 
[0 11] emitting galaxies in the HST COSMOS 2 square de- 
gree field using the Suprime-Cam on the Subaru Telescope 
(Takahashi et al. 2007) as weU as from Ly et al. (2007). The 
overall agreement, both in shape and amplitude, is striking. 
The Mgll rest equivalent width distribution appears to be 
related to the [0 11] luminosity function, and therefore the 
p.d.f. of star formation rate. The connection between Mgll 
absorption and star formation is not only seen as a function 
of redshift, as shown through estimates of the [O n] luminos- 
ity density £011, but appears to also operate as a function 
of Wo. 

Comparing the red curve to the set of blue curves 
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Figure 6. The connection between tlie [On] luminosity function, ^(L), and tlie distribution of Mgll rest equivalent width, dN/dWo- 
Left: ^{L) derived from the distribution of Mgll absorbers (solid red) using the observed scaling relation between [On] emission and 
Mgll absorption, dN/dWo estimated at z ~ 1.2 from the SDSS (Nestor et al. 2005) and a cross-section with a radius of 50 kpc. Three 
observed [On] luminosity functions are shown for comparison. No dust corrections are applied allowing direct comparisons. Right: Mgll 
rest equivalent width distribution at z ~ 1.2 measured with the SDSS (red data points, Nestor et al. 2010). The three blue-colored curves 
show our estimates of d^N/dWodz based on the observed [On] luminosity functions. 



reveals a number of interesting properties: first of all we 
can see that the above relation maps Mgll absorbers with 
1 ^ Wo < 6 A to [O ll] line luminosities ranging from 10*^ 
to about 5 X 10*'^ erg/s, i.e. the range corresponding to the 
exponential decline of the [O ll] luminosity function. Second, 
the match in amplitude shows that the incidence of strong 
Mgll absorbers can be explained by the number density of 
[O ll] bright galaxies and finally the shape of the two dis- 
tributions are in good agreement, further indicating that 
these emission and absorption processes are related. These 
two quantities ultimately allow us to probe the same phe- 
nomenon: star formation. 

In Eq. 12, the terms N /dWo dz, a and £{z) are al- 
ready constrained by direct observations and a set of cosmo- 
logical parameters. Some constraints on the spatial extent of 
the cross-section a and the covering factor C already exist 
but are less robust. Keeping in mind that the scatter in the 
reported values of $(!/) only allows us to compare distribu- 
tions within a factor ~ 2, we now comment on these two 
quantities. Changing the value of the covering factor C in 
Equations 11 and 12 results in shifting the curves vertically. 
By definition, the amount of [O ll] luminosity density traced 
by Mg II absorbers must be lower than the total value, esti- 
mated from direct measurements. Depending on the chosen 
dataset, Figure 6 shows that the amplitude of the [On] lu- 
minosity function implied by Mgll absorbers is either com- 
parable to the measured $(!/) or lower by about a factor 
two. This puts a strong constraint on the covering factor 



Therefore, the covering factor for Mgll absorption around 
star forming galaxies (with Loii 10^^ erg/s) is expected 
to be high, possibly close to unity. Using a different value of 
the cross-section a has several effects: it changes both the 
relation between L and Wo, and the amplitude of the pre- 
dicted shifting the prediction diagonally and changing 
its curvature. We find that only values of 20 < 6max ^ 60 kpc 
do not over-predict the observed luminosity density at cer- 
tain values of Loii- The match between these two sets of 
curves is possible for only a restricted range of impact pa- 
rameters, which encloses the mean value reported by direct 
observations of Mgll-galaxy associations. 

We now look at the relation in Eq. 11 from the opposite 
direction, i.e. we use it as a model for the distribution of 
Mgll rest equivalent width. Using C = 1, as discussed above, 
we can write 



^j^=*(logL)^ 



(14) 



C = C/ Cn =i 0.5 - 1 . 



(13) 



This relation provides a model explaining both the ampli- 
tude, shape and redshift-dependence of the Mg ll rest equiva- 
lent width distribution, with no free parameter. The inferred 
distribution of Mgll rest equivalent widths is shown in the 
right panel of Figure 6 with blue lines using the three mea- 
surements of the [O ll] luminosity function introduced above. 
For comparison, we show the measurement of d?N /dWodz 
from SDSS obtained by Nestor et al. (in prep.) for systems 
centered at z = 1.2, with a half width of 5z = 0.2. This Fig- 
ure allows us to illustrate the accuracy with which the latter 
distribution is measured with current datasets. The estimate 
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of N /dWo dz originates from the analysis of about 45,000 
SDSS quasar spectra spread over a large fraction of the sky. 
It therefore does not suffer from cosmic variance. 

While we have focused on the similarities between the 
shape and amplitude of $(logL) and d^N /AWqAz, we now 
comment on possible departures between these quantities. 
As surveys become larger and selection effects better un- 
derstood, a more robust estimate of the [O II] luminosity 
function will allow us to explore possible deviations from 
the predictions given by Eq. 11. As discussed above, the ra- 
tio between the predicted and observed will provide 
us with an estimate of the covering factor for absorption 
around star forming galaxies. Deviations in shape could also 
indicate variations of the cross-section as a function of [O ii] 
luminosity and/or Mgll rest equivalent width. As mentioned 
above, the overall distribution of Mg ii Wo can be described 
by a Schechter function, with W* ~ 0.6 at z = 1. Observa- 
tional constraints exist down to about 0.1 A. If the corre- 
spondence between Wo and Loii holds in this regime, this 
would provide us with a mean of probing the [O II] luminos- 
ity function down to values significantly fainter than current 
emission derived estimates at these redshifts. 

In summary, we have shown that the distribution of 
Mg II rest equivalent widths and the [O ii] luminosity func- 
tion track each other, as a function of both redshift and 
Wo, suggesting that these two observables probe the same 
phenomenon: star formation. This connection allows us to 
present the first non-parametric model for the distribution 
of Mgii rest equivalent widths, d?N /dWo dz, which appears 
to inherit its amplitude, shape and redshift dependence from 
the [O II] luminosity function or equivalently the p.d.f. of star 
formation in the Universe. 



6 DISCUSSION 

6.1 Absorber-galaxy correlations 

As mentioned earlier, the nature of strong Mgii absorbers 
has been a matter of debate for more than two decades. One 
way to address this question has been to look for observa- 
tional correlations between absorber and galaxy properties. 
However, until recently no compelling detection had been 
reported. This analysis has allowed us to show that Mgii 
absorbers trace a substantial fraction of the overall star for- 
mation rate of the Universe. This property actually sheds 
light on why previous attempts to detect correlations be- 
tween Vyo(Mgll) and galaxy secular properties have either 
failed, or reported rather weat trends (e.g. Kacprzak et al. 
2007). As shown by Brinchmann et al. (2004) at « < 0.2, the 
relation between SFR and galactic stellar mass (measured 
over five decades) is relatively shallow: SFR oc M^". At high 
redshifts, observational constraints on the SFR-stellar mass 
relation are limited to massive galaxies (Noeske et al. 2007, 
and Daddi et al., 2007) but indicate a large scatter in SFR 
at a given mass. Therefore selecting galaxies with respect 
to their SFR, or equivalently Wo, will result in a sample 
with a wide range in stellar mass, broad band luminosity 
and colour. Correlations between absorber rest equivalent 
widths and galaxy secular parameters such as broad band 



luminosity or mass are thus expected to be weak, in agree- 
ments with results from Steidel et al. (1997); Zibetti et al. 
(2007); Menard & Chelouche (2009); Tinker & Chen (2008) 
and Gauthier, Chen, & Tinker (2009). 



6.2 The nature of Mgii absorbers 

The nature of Mg ii absorbers has been long been debated 
and so far no consensus has been reached regarding the phys- 
ical process governing the observed equivalent widths. In this 
work we have shown the connection between the presence of 
Mgll absorbers in galactic halos and neaxby stax formation. 
Gas around a galaxy can be related to star formation in two 
different ways: it can be infalling and feeding an episode of 
star formation or it can be outflowing, due to star forma- 
tion feedback processes such as supernova winds or radiation 
pressure. 

Enhancements of star formation are expected to be trig- 
gered on a timescale of order the dynamical time (t^ = 

R/vc ~ 10**yr) following gas accretion. The time required 
for these systems to travel a distance r of 20 kpc is 

M = mMyx(-^\( j\ , (15) 

^ V20kpcy V300kms-V 

where the value of the characteristic velocity is motivated 
by observations of blue-shifted Mg ii absorption in the spec- 
tra of star-forming galaxies (Weiner et al. 2009) or using 
the observed Mgll rest equivalent width as a proxy for the 
projected velocity: Au ~ 120 (Wo/A) km/s (see section 2.1). 
Such velocities are however comparable to the circular ve- 
locity of L* galaxies and can also be reached by infalling 
material. The above relation indicates that the timescale for 
the gas to travel is comparable to the timescale of a star 
formation episode. Note that we would not expect any cor- 
relation between Wo and instantaneous star formation rate 
if Td < At. 

The high covering factor around star forming galaxies 
implied in our analysis allows us to shod light on the di- 
rection of motion of the gas. Although there is currently 
no direct observational evidence for infalling cold gas onto 
galaxies, results from numerical hydrodynamic simulations 
suggest that gas accretion leading to efficient star forma- 
tion mostly occurs through cold flows, i.e. dense filamentary 
structures penetrating through the diffuse gaseous halos of 
galaxies (Keres et al. 2005; Dekel et al. 2009). Such models 
suggest that these elongated structures give rise to covering 
factors of order 10-20%. In Section 5 we showed that, in or- 
der not to over-predict the observed [O II] luminosity func- 
tion, the covering factor for Mgll absorption around star 
forming galaxies is necessarily high, with C = Cf Cn > 0.5. 
While such a value appears somehow too high for cold flows, 
it is naturally reproduced by outflows, which are known to 
take place with a large opening angle (Veilleux, Cecil, Bland- 
Hawthorn 2005). 

Another argument in favor of outflows comes from the 
similarities between the results presented in this study, i.e. 
the connection between intervening Mg ii absorbers and star 
formation, and the recent results by Weiner et al. (2009) 
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who showed that blue-shifted Mg ii absorption, i.e. outflow- 
ing gas, is ubiquitous in star- forming galaxies at a: ~ 1. 
These authors found the presence of outflows in all subdi- 
visions of their sample of star-forming galaxies, covering a 
factor ten in SFR and B-band luminosity, and a factor 30 
in stellar mass. Their blue-shifted Mgii systems appear to 
share all the properties of the intervening systems used in 
our study: rest equivalent widths, redshift range, high cov- 
ering factor and especially the correlation with star forma- 
tion. Occam's razor (c. 1320), which tells us that entities 
should not be multiplied unnecessary, leads us to postulate 
that the nature of the Mg ll gas is the same in both types of 
analyses. In other words, blue-shifted Mg ll absorbers seen in 
self absorption and intervening Mg ll absorbers might be the 
same gas clouds. This implies that the blue-shifted absorp- 
tion seen in star forming galaxies (e.g. Heckman, Armus, & 
Miley 1990; Martin 1999; Weiner et al. 2009 may well extend 
to radii up to tens of kpc. 



6.3 A new probe of star formation 

Our analysis has shown that the bulk of star formation is 
traced by Mgii absorbers. This connection provides us with 
a new tool to explore star formation in the Universe, over 
cosmological times. Here we discuss its potential and high- 
light differences with existing techniques. 

• Absorption measurements offer a number of advan- 
tages: the detectability of absorption lines does not depend 
on redshift (while galaxy fluxes drop rapidly at high red- 
shifts) and perhaps more importantly, the observed proper- 
ties of absorption lines are not sensitive to dust extinction. 
Mg II absorbers can therefore provide us with a powerful and 
complementary tracer of star formation. 

• Absorber systems allow us to apply emission line studies 
in a noise-dominated regime. Only a handful of low-redshift 
Mg II absorbers allow for the direct detection of [O ll] emis- 
sion lines in individual systems. For the vast majority, no 
direct detection is possible in SDSS quasar spectra. How- 
ever, being able to average over thousands of lines-of-sights 
allows us to estimate the mean [O ll] luminosity of the corre- 
sponding population, down to emission levels comparable to 
what is currently obtained with some of the deepest narrow 
band surveys. 

• This technique does not suffer from contamination from 
misidentified emission lines, a severe limitation for narrow- 
band imaging surveys of [O II] emission (see Takahashi et al. 
2007; Ly et al. 2007). 

• In optical spectra Mgii absorbers are detected from 
the ground over the range 0.3 < z < 2.2. Irffrared spec- 
troscopy can extend this redshift path up to the reioniza- 
tion epoch. For example, the upcoming near-infrared spec- 
trograph FIRE (Simcoe et al. 2008) will allow the detec- 
tion of Mg n absorption and [O ll] emission up to a ~ 6, 
for which the rmmber of quasars will become the limiting 
factor. Infrared spectroscopy will allow us to map the star 
formation history beyond redshift two, with the same tracer, 
and provide a new set of constraints on the decline of the 
SFR density at high redshift. 



A vast amount of Mg ii data has been accumulated over 
the years. Its analysis has mainly focused on attempting to 
reveal the nature of these systems. We can now use these 
datasets to obtain new insights into star formation, over 
cosmological times. About 16,000 Mg ii absorbers have been 
detected in SDSS quasar spectra Quider et al. (2010). The 
analysis of the entire sample (DR7) will double this number. 
The details of Mg ii absorbers (number of components, ve- 
locity widths, etc.) accessible through high-resolution spec- 
troscopy can be used to put constraints on star formation 
processes. Understanding this type of information might al- 
low us to constrain the recent star formation history of 
galaxies. Such a technique will be complementary to spectral 
energy distribution analyses. 

More work needs to be done in order to understand 
the origin of the empirical [Oll]-Mgll scaling relation and 
especially its redshift dependence. If the redshift depen- 
dence of ELon(W^o) can be explained, Mgll absorbers can 
then, alone, be used to trace star formation in the Uni- 
verse. Finally, comparing SFR(2;) from emission and absorp- 
tion techniques might shed light on the redshift dependence 
of the dust extinction affecting emission based estimators. 
With appropriate wavelength coverages, similar absorption- 
emission scaling relations can be explore, for example be- 
tween Mgll absorption and H-a emission. It will be inter- 
esting to extend the current analysis to other absorption 
transitions (Fell, CIV, etc.) and test their connection to the 
process of star formation. 



7 SUMMARY 

We present an empirical connection between cold gas in 
galactic halos and star formation. Using a sample of more 
than 8,500 Mgll absorbers from SDSS quasar spectra, we 
report the detection of a 15 a correlation between the rest 
equivalent width Wo of Mg ii absorbers and the associated 
[O ii] luminosity, an estimator of star formation rate. This 
correlation allows us to show that 

(i) Mg II absorbers trace a substantial fraction of the [O ll] 
luminosity density in the Universe, i.e. global star formation. 
The SFR density probed by Mgll absorbers recovers the 
overall star formation history of the Universe derived from 
classical emission estimators up to « ~ 2. 

(ii) the distribution function of Mgll rest equivalent 
widths, dN/dWo inherits both its shape and amplitude from 
the [On] luminosity function $(!/). These distributions can 
be naturally connected, without any free parameter. 

(iii) the covering factor for cold gas around star forming 
galaxies appears to be high: C > 0.5. The determination of 
this value is currently limited by the scatter in the reported 
estimates of [O n] luminosity densities. 

These results not only shed light on the nature of Mg ll ab- 
sorbers but also provide us with a new probe of star forma- 
tion, in absorption, i.e. in a way which does not suffer from 
dust extinction and with a redshift-independent sensitivity. 
As shown in this analysis, such a tool can be applied in a 
noise-dominated regime, i.e. using a dataset for which emis- 



© 0000 RAS, MNRAS 000, 000-000 



12 B. Menard et al. 



sion lines are not detected in individual objects. This is of 
particular interest for high redshift studies. 

We argue that outflows appear as a favored explanation 
regarding the nature of Mgii absorbers and that blue-shifted 

Mgll absorption seen in the spectra of star forming galax- 
ies could be the same systems, implying that the observed 
outflowing gas can reach radii of ~ 50 kpc. 
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